Trapping Dynamics with Gated Traps: Stochastic Resonance-Like Phenomenon 
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We present a simple one-dimensional trapping model prompted by the problem of ion current 
across biological membranes. The trap is modeled mimicking the ionic channel membrane behaviour. 
Such voltage-sensitive channels are open or closed depending on the value taken by a potential. Here 
we have assumed that the external potential has two contributions: a determinist periodic and a 
stochastic one. Our model shows a resonant-like maximum when we plot the amplitude of the 
oscillations in the absorption current vs. noise intensity. The model was solved both numerically 
and using an analytic approximation and was found to be in good accord with numerical simulations. 



During the last decade we have witnessed a growing interest in the phenomenon of stochastic resonance (SR). 
This phenomenon is related to the enhancement of the response of a nonlinear system due to the interplay between 
nonlinear oscillations and noise. SR has been studied in different physical, chemical and biological contexts [EJ. In 
particular, it has been found to play a relevant role in several problems in biology: mammalian sensory systems, 
increment of the tactile capacity, visual perception, low frequency effects and low amplitude electromagnetic fields, 
etc@. 

Among the many experiments showing SR there is one related to the measurement of the current through voltage- 
sensitive ion channels in a cell membrane ||. These channels switch randomly between open and closed states, 
controlling the ion current. This and other related phenomena have motivated several theoretical studies of the 
problem of ionic transport through biomembranes. Those studies have used different approaches, as well as different 
forms of characterizing stochastic resonance in such systems Q . 

In this work we want to analyze the effect of the simultaneous action of a deterministic and a stochastic external 
field on the trapping rate of a gated imperfect trap. Here we do not pretend to make a precise modeling of the 
behaviour of an ionic channel. We propose a simple model of dynamic trapping. This model is motivated by the 
above indicated experiment on cell membranes |3j. Our main result is that even a simple model of a gated trapping 
process can show a SR-like behaviour. 

The study of gated trapping processes, i.e. a trapping process where the traps have some kind of internal dynamic 
has attracted considerable interest || . Many authors discussed the way to link the gated trapping processes with the 
measured behavior of the so called ionic pumps Q. For example, among other factors, the ion transport depends on 
membrane electric potential (which plays the role of the barrier height) and can be stimulated by both dc and ac 
external fields. 

We based our study on the so called stochastic model for reactions , that has been generalized in order to include 
the internal dynamics of traps. The dynamical process consists of the opening or closing of the traps according to 
an external field. Such a field has two contributions, one periodic with a small amplitude, and the other stochastic 
whose intensity will be the tuning parameter. 

Our starting model equation is 

d t p(x, t) = Ddlp{x, t) - j(t)S(x)p(x, t) + n u , (1) 

where 7 is a stochastic process that represents the absorption probability of an individual trap in our system, p is the 
particle density (particles that have not been yet trapped); for a given realization of 7, x is the coordinate over the 
one-dimensional system and n u is a source term that represents a constant flux of ions. The injection of ions can be 
at a trap position or at any other position. In this last case the ion can diffuse to the trap position. This diffusion 
coefficient would represent an effective diffusion through the volume rather a diffusion over the membrane surface. 
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Each 7 is modelled as follows: 



7(i)=7*0[SBin(wt) + 



(2) 



where 0(x), the Heaviside function, determines when the trap is open or closed. The trap works in the following 
manner: if the signal composed of the harmonic part plus £ (the noise contribution) reaches a threshold £ c the trap 
opens, otherwise it is closed. We are interested in the case where £ c > B, that is, without noise the trap is always 
closed. When the trap is open the particles are trapped with a given frequency (probability per unit time) 7*. In 
other words the open trap is represented by an "imperfect trap" . Finally, in order to complete the model, we must 
give the statistical properties of the noise £. We assume that £ is an uncorrelated Gaussian noise of width £ , i- c - 



St < 2 



(3) 



Note that this is not a standard white noise [|| due to the Kronecker symbol (instead of Dirac's delta 5(t — u)). 
For future reference we compute {"f(t)} (brackets mean averages over all realizations of the noise). The result is 



-erfc 



(£ c -5sin(wi)) 



(4) 



Setting n(x, t) = (p(x,t)) and averaging Eq. (nh we can write 



n(x, t) = n u t 



duJ(u)G(x, t — u), 



(5) 



where G(x,t) = Go(x,t) is the free diffusion propagator, and J(t) 
This current satisfies the relation 



(lj(t)p(jl,t)) is the current trough the trap. 



J(t) = <7(*)> 



n u t — / duJ(u)G(0,t — u) 




(6) 



We solve this equation numerically, using a standard "trapezoidal" integral approach. 

The present problem can be extended to a finite lattice with periodic boundary conditions. In this case Eq. (|l|) is 
replaced by 



d t p(x, t) = Dd 2 x p{x, t) 



E 



~f(t)5(x - jl) + n u , 



(7) 



where, due to the boundary conditions, the sum runs over all the images of the trap. The current satisfies the same 



Eq. 



but now G(x, t) — J7 Gq(x — jl, t). This equation can be solved numerically. 



In addition to the numerical solution, as in typical SR studies [[TJ , we can make some analytic calculations using a 
perturbation procedure. We assume B/£q to be small and expand (j(t)) (or {~fj{t)}) in the following form 



(7(*)> = (7i(*)> =7o+7i sin(wt) + 0((B/Ho) 2 ), 



(8) 



where 



7i 



—erfc 


( ?C 1 


2 




7*£ 


( ? 




\/2<o 



(9) 
(10) 



Now we can transform Eq. (||) to the Laplace domain and solve it iteratively. The expressions obtained are rather 
complicated and for simplicity we present only the asymptotic behavior. For an infinite lattice the asymptotic current 
is 



J(t) — 4n„ 



— + \C\ cos{u;t - 

7T 



(11) 



where 
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C= \C\exp(i</>) = —j=— r- , (12) 

7o(7o/(2V-D) + V*w) 

while for a finite one it is 

J(t) ~ nj + \C\ cos(ujt + </>), (13) 

where 



C = \C\exp{i<j>) = — j— T — . (14) 



7oV«w< 



The qualitative behaviour of the system changes since for an infinite lattice the current grows without limit while in 
a finite lattice the current reaches a stationary value when the source particle supply is balanced with the trapping 
rate. 

These expressions are valid only for large values of the width noise £o- In the opposite case, i. e. for small values 
of £o> the expansion of (7) (or (79)) as given in Eq. (ph is not appropriate. In this small £0 limit we note that J(t) is 
small and then we can obtain the following approximate solution for Eq. (^J) 

J(t) ~ (l{t))n u t, (15) 

that is independent of G and hence is valid for both situations (finite and infinite lattice) . This expression is valid for 
short times such that we can neglect the integral in Eq. (^). 

We choose to quantify the SR-like phenomenon by computing the amplitude of the oscillating part of the absorption 
current given by: 

= J\mn(ut)=l ~ ^lsin(u>t)=-l (16) 

The previous limits allow us to explain the qualitative behaviour of the system as follows. For a small noise intensity 
the current is low (remember that £ c > B), hence A J is small too. For a large noise intensity the deterministic part 
of the signal (the harmonic one) becomes irrelevant and the AJ is also small. Therefore, there must be a maximum 
at an intermediate value of the noise. 

We note that the case corresponding to a periodic array of fixed traps only requires replacing 7 for jj in the sum 
of Eq. (0). In the case of the jj being independent of each other, i.e. Jj(t) = 7*#[i? sin(w£) + £j — £ c ], where 
(£j(t)£k{u)) = Sj,kSt.u^o the current within the SM is exactly the same that for the finite lattice case. 

The simulations were performed on one dimensional lattice of L sites with periodic boundary conditions. Initially 
there are no particles on the lattice. The particles are injected randomly every l/(n u L) units of time, with uniform 
distribution over the lattice, and are allowed to perform a continuous time random walk (characterized by the jump 
frequency at each neighbor site q). There is no restriction on the number of particles at each site. A particle can be 
removed from the system with a given probability distribution characterized by 7 when it reaches the trap site. A 
detailed description of the algorithm used can be found in Q . The reaction times were generated according to the 
following probability density function 

p(t) = exp(- /*<7(0)dt , V (17) 







All simulations shown in the figures correspond to averages over 1000 realizations. 

The results of both, numerical evolution and simulations, are shown in Figs. 1 to 3. In Fig. 1 we show the 
temporal behaviour of the absorption current as a function of time. Figure la shows the evolution obtained from 
numerical integrations for a single trap (infinite system) and a periodic array of traps in a time interval such that the 
periodic array reaches its stationary behavior ||. It is worth remarking that a system with a single trap does not 
have a stationary behaviour. However, the figure can induce the false impression that the single trap also presents 
a stationary behaviour. This is only an artifact of the short time period shown. The insert of Fig. la shows, in a 
much longer time scale, such a non stationary behaviour for the single trap problem. The agreement between both 
theoretical results makes it difficult to distinguish among them, particularly at short times. Figure 16 depicts the 
results obtained from the simulations (points) and the corresponding numerical integration of the single trap (solid 
line) for a fixed value of £o- One see the good agreement between the two curves. When looking at such curves, and 
for different values of the noise intensity, it is apparent that for intensities smaller or larger than an optimal one, the 
amplitude is reduced, while in its neighborhood the amplitude becomes a maximum. This is better seen when plotting 



3 



A J, the amplitude of the oscillating part of the absorption current as a function of the noise intensity, as can be seen 
in Fig. 2 for two different times. Such results make apparent the existence of a maximum as a function of £o- Here 
we find good qualitative agreement between the results obtained by numerical integration of Eq. (^) and those from 
simulations. The insert of Fig. 2 shows the results for the initial slope indicated in the figure by h, corresponding 
to the kind of curves shown in the principal part of the figure, as a function of time. That is, the slope of AJ for 
low noise intensities and short times. Both results, obtained by numerical integration and simulations, show linear 
dependence with time as predicted in Eq. (15), however the actual values of the slopes are sightly different. 

Finally in Fig. 3 we depict the result for <j>, the phase shift between the input periodic signal and the absorption 
current (as indicated in Eqs. (12) and (14)). We have compared the results from the numerical integration of Eq. (||) 
and those from simulations, finding good agreement. 

The present results show that even a simple trapping process with an adequate internal dynamic can present a SR 
like phenomenon. This suggests the possibility that such a kind of phenomenon could be ubiquitous in more complex 
trapping like processes. 

A more elaborate model aimed at describing the experiment in Ref. || , including a more realistic representation of 
the ion channels is under study and will be the subject of further work. 
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FIG. 1. Temporal evolution of the absorption current of particles for a gated trapping system, a) Time evolution obtained 
from numerical integration for a single trap (solid line) and for a finite system with periodic boundary conditions (points). We 
show the stationary behaviour for the periodic array of traps. The non stationary behaviour of the single trap problem, as 
discussed in the text, is apparent from the insert, b) Time evolution obtained from simulations (points), and from numerical 
integration of a single trap system (solid line). The parameters are: u = 0.1, n u — 0.01, £ c = 1.2, B — 1.0, qA = 0.1, 7* = 10, 
£0 = 0.3 and L = 100. 



FIG. 2. Variation of the final value of AJ (amplitude of the oscillating part of the absorption current) as a function of £o- 
The results obtained from the numerical integration of Eq. (6) are indicated by solid lines, while those from simulations by 
circles. We have chosen the following time parameters; a)t — 517 and b)t = 895. The insert shows the initial slope h as a 
function of time. The results obtained from the numerical integration are indicated by triangles and those from simulations 
by circles. The linear fitting is indicated by continuous and dashed lines respectively. The results were obtained for low noise 
intensity values (£0 < 0.1). 
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FIG. 3. Phase shift as a function of £o- The results obtained for the numerical integration are indicated by a solid line 
while those from simulations by circles. 
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